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Intracellular transport based on molecular motors and its regulation are crucial to the functioning 
of cells. Filamentary tracks of the cells are abundantly decorated with non-motile microtubule- 
associated proteins, such as tau. Motivated by experiments on kinesin-tau interactions [Dixit et 
al. Science 319, 1086 (2008)] we developed a stochastic model of interacting single-headed mo- 
tor proteins KIF1A that also takes into account the interactions between motor proteins and tau 
molecules. Our model reproduce experimental observations and predicts significant effects of tau on 
bound time and run length which suggest an important role of tau in regulation of kinesin-based 
transport. 
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I. INTRODUCTION 

Intracellular transport is fundamental for cellular func- 
tion, survival, and morphogenesis. Kinesin superfam- 
ily proteins (also known as KIFs) are important molec- 
ular motors that directionally transport cargoes along 
microtubules (MTs), including membranous organelles, 
protein complexes and messenger ribonucleic acids (mR- 
NAs) [3| . Disruptions or defects of MT- based trans- 
port are observed in many neurodegenerative diseases 
[3, Q . MTs are decorated with non-motile microtubule- 
associated proteins (MAPs) that promote MT assembly 
and play important roles in organizing the MT cytoskele- 
ton [J, Q. Kinesin proteins interfere with MAPs, being 
the latest able to inhibit active transport of cytoplasmic 
material. Kinesin competes with MAPs to bind to the 
MT surface and MAPs bound to MTs might also block 
the path of motor proteins. Tau is a mainly neuronal 
MAP, enriched in the axonal compartment [6| that has 
been shown to inhibit plus end-directed transport of vesi- 
cles along MTs by kinesin Q. Tau reduces not only the 
attachment frequency of kinesin to MTs but also the dis- 
tance that kinesin travels along the MT in a single run 
Q. It is also known that when single kinesin motors 
encounter tau patches on the MT, most of the motors 
detach from the MT surface Q. 

Some recent theoretical models of interacting molec- 
ular motors [lol - fl3| are extensions of asymmetric simple 
exclusion processes (ASEP) [TJ] in which motors are rep- 
resented by particles that hop along a one-dimensional 
lattice with hard-core exclusion. A paper by Parmeg- 
giani et al. [H[ introduced an ASEP-like model that 
relaxes the restriction of the conservation of particles in 
the bulk, allowing the attachment or detachment of par- 
ticles in any site. In this paper wc consider specifically 
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the effect of MAP tau in the intracellular traffic of single- 
headed kinesin motor protein KIF1A [lj|- A model for 
KIF1A motor transport has been proposed by Nishinari 
et al. [TtJ that enriched Parmeggiani's model via two sig- 
nificant modifications. Firstly, the model explicitly incor- 
porated the Brownian ratchet mechanism for individual 
KIF1A motors, i.e., the mechanochemical coupling in- 
volved in their directed motion, including the adenosine 
triphosphate (ATP) hydrolysis that fuels the cyclic steps 
of the motor. Secondly, the model included a construc- 
tion based on parameters that have direct correspondence 
with experimentally controllable quantities, of which the 
two most important are ATP and motor concentrations. 
In contrast to earlier models of molecular motors intra- 
cellular transport fIol - fl3l [Tol [l7| . which consider only 
motor-motor interactions, our model incorporates the ef- 
fect of tau on kinesin motors dynamics. One remarkable 
exception is the paper by Chai et al. |l8| which con- 
siders transport of molecular motors in the presence of 
static defects. Unlike the model investigated by Chai et 
al., we consider specifically MAP tau effects in kinesin 
dynamics via parameters that are related to experimen- 
tally observed quantities by Dixit et al @. 



II. LATTICE-GAS MODEL FOR KIF1A 
DYNAMICS IN PRESENCE OF TAU 
MOLECULES 

We work with a onc-dimcnsional lattice of L sites that 
represents one MT's protofilament. Each site on the lat- 
tice corresponds to one motor-binding site on the protofil- 
ament. The lattice parameter is taken to be 8nm, which 
is the separation between adjacent binding sites of the 
MT. Each KIF1A is represented by a particle that, when 
bound to one site of the filament, can hop to any of the 
two nearest-neighbor sites. The motor protein can also 
attach to (detach from) the MT; this is modeled as a par- 
ticle creation (annihilation) on the filament. Each site of 
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the filament has two internal states that represent the 
two possible bound states in which a KIF1A motor can 
be related to the MT. In the state 1, the motor protein 
is strongly bound to the MT, whereas in state 2 it is 
weakly bound (the motor protein is tethered to the MT 
by an electrostatic attraction that prevents it from dif- 
fusing away from the filament) flrjj . Therefore, there are 
three possible states for each site: if the site is empty, 
1 if the site is occupied with a motor in state 1, and 2 if 
the site is occupied with a motor in state 2. We consider 
a distribution of tau molecules decorating the filament. 
Each binding site of the filament can either have one 
tau molecule or no molecule at all. To characterize the 
amount of tau present in the MT, we will refer to the 
concentration of tau molecules, denoted by Ctau, defined 
as the number of tau molecules per site in the filament. 
The tau concentration therefore ranges between and 1, 
with these extreme cases representing no tau molecules 
attached to any site of the MT and a tau molecule at- 
tached to every site of the filament, respectively. 

The dynamical evolution of the system is described 
with transition rates that reflect the stochastic nature of 
the movement of the motor protein and of its interaction 
with tau molecules. We run Monte Carlo simulations 
that update the state of all sites of the filament following 
the random sequential method. In one Monte Carlo Step 
(MCS) we give to all sites, in a random order, the oppor- 
tunity to change their state of occupation. The change 
of state of a given site is regulated by transition rates 
which depend on the state of occupation of the site, the 
state of occupation of its nearest-neighbor sites, and also 
depend on whether or not there are tau molecules in the 
sites under consideration. 

Figure[T]shows a scheme with the updating possibilities 
for a site in the bulk (different rules apply for the two 
ends of the lattice and they are specified below). Let 
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FIG. 1. Scheme of the model showing the updating possibil- 
ities and the associated transition rates for a bulk site. The 
effect of tau, preventing a motor attachment or forcing a mov- 
ing motor to detach is graphically depicted as branching ar- 
rows. 



that a tau molecule prevents the attachment of a KIF1A 
motor protein to site i. It is equal to if there is no 
tau molecule in the site, and equal to p a if there is a 
tau molecule in the site. If the site i is occupied with 
a motor in state 1, two things can happen. Firstly, the 
motor protein can detach from the site, with a rate Ud, 
and secondly, the motor protein can hydrolyze the ATP 
molecule and consequently change to state 2, with a rate 
uih- If the site is occupied with a motor in state 2, several 
possibilities may arise. The motor protein can diffuse 
to one of its nearest-neighbor sites with a rate uj (1 — 
T m {i — 1)) towards the minus end, and u>b(l — T m (i + 1)) 
towards the plus end (staying in state 2); it can release 
adenosine diphosphate (ADP) and move forward with a 
rate w/(l — T m (i)) (binding to the next site in state 1), 
or, finally, it can release ADP and stay in the same site 
with a rate oj s (changing to state 1). T m (i) represents the 
probability that a tau molecule in site i forces a moving 
motor protein to detach. It is equal to if there is no 
tau molecule in the site, and equal to p m if there is a tau 
molecule in the site. In this updating process a change 
of site is possible only if the target site is empty. 

For the dynamics at the ends, we take a and 5 instead 
of uj a as the attachment rates for the minus and plus 
ends, respectively. In the same way, 71 and (3\ (instead 
of uid) are the detachment rates, and 72 and ($2 (instead 
of u)b) are the exit rates of the motors due to Brownian 
motion for the minus and plus ends, respectively. 

Rate uj means that in an infinitesimal time interval dt, 
the probability of the event occurring is uidt. 



A. Mean-field equations 

Let us denote by r, and qi the probabilities of finding a 
KIF1A motor in states 1 and 2, respectively, at the lattice 
site i at time t. The master equations for the dynamics 
of motors in a bulk site are given by: 

dv ' 

-7T = (1 - T a (i)) UI a (1 - Ti - qi) - (uj d + Uh)ri + UlsQi 
at 



+ (1 - r m (i)) wjqi-i (1 - n - q t ) , 

dqi n . 

— = uj h r t - uj s q t - tu f q l (1 - r t+x - q i+l ) 



(1) 



- CJ b qi (2 - n-i - q t -i - r. l+1 - q i+1 ) 

+ (1 - r m (£)) u b (gj_i + q l+ i) (1 - n - qi) , (2) 



where 



and 



p a if there is a tau molecule on site i 
if there is not a tau molecule on site i ' 

(3) 



us consider the ith site: if the site is empty, the only 
possible state change is that a motor binds to it in state 1, 
with a rate w a (l — T a (i)). r a (i) represents the probability 
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p m if there is a tau molecule on site i 
if there is not a tau molecule on site i 



(4) 
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The corresponding equations for the minus and plus 
ends can be written in a similar way, with their appro- 
priate constants as described previously. 



B. High-density limit 

We can obtain an approximate solution in the high 
oj a limit, which corresponds to a high motor concentra- 
tion regime with a jammed region that covers the whole 
MT [191 ]. Assuming periodic boundary conditions, the 
solutions (ri,qi) = {r,q) of the mean- field equations in 
the steady state, with the additionally assumption that 
T a (i) = r a Ctau and r m (i) = r m Ctau, are found to be 



l + (K s (l-n f )-l)g 
l + K d + K h {l-n f ) ' 



= -A - y/A 2 - 414 Wd + Wh + K s )(l - Sl f )j 
q 2(K d + (K h + K s )(l-Q f )) 

where A = K h {Q, } - 1) - K d (tt s + 1) - (f2 s + h ), 
K — ^ Ki — ^ Kj 



(5) 



(6) 



(l-Ta)w. 



O 



(l-r a )w 



(l-T a )Wa ' 



A system with periodic boundary conditions is not re- 
alistic. However it provides a good approximation in the 
high motor concentration limit with open boundary con- 
ditions. Figure [2] shows a comparison between these an- 
alytical expressions and numerical simulations for three 
cases corresponding to different magnitudes of the tau- 
kincsin interactions. 

The limit with no tau molecules present in the cytosol 
can be obtained taking r a — and r rn = 0. In that 
case, expressions §5§ and ^ reduce to the ones found 
by Nishinari et al. in the model for KIF1A dynamics 
without tau molecules [l7j |. 



III. BOUND TIMES AND RUN LENGTHS 

We work with a MCS representing 0.71ms of real 
time. Every simulation has a fixed concentration of tau 
molecules and starts with no motor attached to the fil- 
ament. In order to construct the distribution of tau 
molecules with the desired concentration, we give to each 
site of the filament the chance of having a tau molecule 
with probability Ct au . Dixit et al. found that tau patches 
were stable over the time course of several minutes while 
the average duration of motor protein runs are in the or- 
der of seconds or tens of seconds [j| . Therefore, the tau 
distribution is held constant during the complete course 
of the simulation. Our model captures the essence of the 
experimental results on tau-kincsin interactions found by 
Dixit et al. @, that can be summarized as: i) the fre- 
quency of kincsin binding to MT depends on local tau 



concentration; ii) the average binding frequency of ki- 
nesin in presence of tau in the cytosol reduces to approx- 
imately one third of the value without tau; iii) approx- 
imately half of the times that a kinesin motor protein 
encounters a tau molecule, it detaches from the MT, the 
remaining half times it passes or pauses. Therefore, the 
probability that a tau molecule prevents a kinesin motor 
protein binding to the MT is taken to be p a = 0.67, and 
the probability that a tau molecule forces a moving ki- 
nesin motor protein to detach from the MT is chosen to 
be p m = 0.5. All parameters that have no relation with 
tau were chosen following Ref. (l7j . 

For the cases shown in Figs. [3] and [4] the model pa- 
rameters were taken to be L = 600, p a = 0.67, p m = 0.5, 
io a = a = 0.00001 [ms- 1 ], uj d = ft = /3 2 = O.OOOlfms" 1 ], 
uj s = 0.1451ms- 1 ], uj = 0.0551ms- 1 ], u b = O^ms" 1 ], 
and 7i = 72 = S = 0. The observables shown were ob- 
tained averaging the results of 1000 simulations with the 
same parameters except for the random number genera- 
tor seed. The error associated to the observables corre- 
spond to the standard deviation. Each one of the simula- 
tion generates a new tau distribution (with the same con- 
centration) and evolves for 2xl0 6 MCSs after a thermal- 
ization of 2xl0 5 MCSs which guarantees that the system 
reaches the stationary state. Figure |3] shows the behavior 
of the bound time as a function of uj^. In each simula- 
tion the bound time is calculated by averaging the inter- 
vals between attachment and detachment of each KIF1A. 
It is worth remembering that u>h is a monotonous in- 
creasing function of ATP concentration (ojh = 0.25ms _1 
corresponds to saturating ATP concentration) [T3|. For 
Ctau = there are no tau molecules present in the cy- 
tosol (in this case, the values of p a and p m are irrele- 
vant) and the bound time grows with u>h from near 15s 
for low ATP concentration to almost 25s for saturating 
ATP concentration. When there are no tau molecules 
bound to the MT, motor proteins can only detach from 
the filament when they are in state 1. So, as oj% in- 
creases, the ATP concentration increases and the time 
that a motor is in state 2 becomes larger (if there are 
enough ATP molecules, the motor hydrolyzes one and 
changes to state 2 almost immediately). Thereby, the 
time duration of the motor runs in the MT increases. On 
the other hand, in the presence of tau molecules, KIF1A 
motors can be forced to detach even if they are in state 
2. When a moving motor protein (either actively mov- 
ing forward or diffusively moving towards either end) en- 
counters a tau molecule in its target site, it is forced to 
detach from MT with a probability p m . In this way, the 
presence of tau molecules shortens the averaged bound 
time. The greater the values of tau concentration, Ctau-, 
the larger the number of KIFlA-tau encounters and the 
stronger the effect of tau on the averaged bound time. 
This can be observed in Fig. [31 where the values of av- 
eraged bound time are ordered with increasing value of 
Ctau from top to bottom. The results for C ta u = 0.005, 
in which the tau concentration is very low, are similar to 
the case without tau, i.e., the averaged bound time in- 
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FIG. 2. (Color online) Stationary density profiles for p a = 0.67, p m = 0.5 (left), p a = 0.067, p m = 0.05 (center) , and p a — 
0.0067, p m = 0.005 (right). Model parameters: L = 600, Ctau = 0.01, u; a = a = O-OOlfms" 1 ], iv d = fa = fa = 0.00011ms" 1 ], 
uj s — 0.145[ms _1 ], ujf — 0.055[ms ], = 0.6[ms _1 ], and 71 = 72 = S — 0. The red upper (green lower) points correspond 
to numerical results for the state 1 (state 2) density. The dashed lines are the analytical approximate predictions for periodic 
systems with the same parameters. Note that the chosen uj a value allows us to realize the high density limit 
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FIG. 3. Averaged time between attachment and detach- 
ment of motor proteins on the filament as a function of the 
hydrolysis rate. 



FIG. 4. Average length spanned by KIF1A motors on the 
filament in single runs as a function of the hydrolysis rate. 



creases with uih- For Ctau = 0.01 and Ct au = 0.1 the tau 
concentration is large enough to cancel the effect of the 
high ATP concentration and the averaged bound time 
remains almost constant for all oj% values. Moreover the 
averaged bound time for Ctau =0.1 presents a slightly 
decreasing tendency with increasing values of Uh (aver- 
aged bound time decreases approximately from 0.22s to 
0.13s ) which corresponds to a tau- ruled behavior. A 
high ATP concentration accelerates the dynamic of the 
motors favoring the encounters with tau molecules and 
increasing the chances of tau mediated detachment. 

The effect of tau on KIF1A dynamics considering the 
averaged run length (calculated by averaging the dis- 
tances that each motor travels between attachment and 
detachment) can be established in similar terms. Firstly, 
we can note in Fig. [4] that in absence of tau molecules 
the averaged run length increases with u>h- The ATP 
molecules are the fuel that KIF1A motors use to actively 
move toward the plus end of MT, so that a greater ATP 



concentration accelerates the dynamic of the motors en- 
larging the distance that the motors travel before detach- 
ing from MT. In presence of tau molecules, the movement 
of KIF1 A along the MT can give room to encounters with 
tau molecules so that a competition begins when we in- 
crease the ATP concentration. On the one hand, as we al- 
ready mentioned, the availability of ATP molecules favors 
longer motor run lengths, but, on the other hand, it also 
favors larger numbers of kinesin-tau encounters, which 
shortens motor run lengths. Therefore for Ctau = 0.005 
and Ctau = 0.01 the run lengths are still increasing func- 
tions of tUh but for Ctau = 0.1 it remains almost constant 
for the whole range of ui^. In this case, the tau concen- 
tration is large enough to compensate the higher ATP 
concentration with a greater number of motor detach- 
ments. For the same value of ujh (ATP concentration) 
the run length is always smaller for larger values of Ctau ; 
that is, the effect of tau is to shorten the distance that the 
motors travel along the MT. Figure [5] shows the differ- 
ential effect of several tau concentrations on the KIF1A 
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run length distributions. Note that these results agree 
well with the experimental observations communicated 
by Dixit et al. in Rcf. Q. 
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FIG. 5. Histograms that show the normalized number of 
events for the run length of KIF1A motors. Model pa- 
rameters: L = 600, p a = 0.67, pm = 0.5, tvd = Pi = 
/3 2 = 0.00011ms- 1 ], uj s = 0.145[ms" 1 ],Lj / = 0.0551ms- 1 ], 
tub — 0.6[ms~ ], 71 = 72 = 5 = 0, ojh = 0.2[ms ] and 
uj a = a = O.OOOOltms" 1 ]. Total number of MCSs = 2xl0 6 , 
after a thermalization of 2x10° MCSs. 



IV. CONCLUSIONS 



tend to preserve one of its virtues, that is, to relate the 
properties of the transport to experimentally controllable 
quantities. In this way, we model the effect of one tau 
molecule by means of two parameters p a and p m . The 
effect the tau molecule has in preventing kinesin motor 
binding to MT is characterized by p a , and the effect it 
has forcing moving motor proteins to detach from the 
MT upon an encounter, by p m . These parameters values 
have been determined from experimental results obtained 
by Dixit et al. @ and have proven to be robust from the 
point of view of the simulation. That is, small changes on 
the values of these parameters do not modify the quali- 
tative behavior of the observables discussed here. More- 
over, changes in orders of magnitude on the parameters 
p a and Pm-, change the value of the observables in the 
same order of magnitude (result not shown). Our model 
coincides with the model by Nishinari et al. [l?} in the 
limit where no tau molecules are present in the filament, 
and the same is true for the analytical approximate high 
density results shown in Eqs. [S]and|ni The processivity 
of molecular motors, which is crucial for the intracellu- 
lar transport, can be defined in several ways [20| . two 
of that being the attachment lifetime of the motor pro- 
teins in the filament (bound time) and the mean length 
traveled by the motor proteins in a single run in the fila- 
ment (run length). We found significant effects on both 
bound time and run length of KIF1A motors due to the 
presence of bound tau molecules in the MT. The results 
shown here suggest that tau could play an important role 
in regulation of intracellular transport by kinesins. 



In this paper we have presented a novel stochastic 
model for the intracellular transport by KIF1A motor 
proteins that takes into account the effect of the MAP 
tau. Our model is an extension of the model of Nishinari 
et al. for the trafficking of KIFIA without tau, and in- 



ACKNOWLEDGMENTS 

J.S. and P.W.L. want to thank SECyT-UNC (Ar- 
gentina) for financial assistance. 



[1] J. Howard, Mechanics of Motor Proteins and the Cy- 

toskeleton (Sinauer Associates, Sunderland, MA, 2001) 
[2] M. Aridor and L. A. Hannan, Traffic 1, 836 (2000) 
[3] I. Toyoshima, K. Kato, M. Sugawara, C. Wada, and 

O. Masamune, Neuroscience Letters 249, 103 (1998) 
[4] G. Lee, Current Opinion in Cell Biology 5, 88 (1993) 
[5] N. Hirokawa, Current Opinion in Cell Biology 6, 74 
(1994) 

[6] S. C. Papasozomenos and L. I. Binder, Cell Motility and 
the Cytoskeleton 8, 210 (1987) 

[7] A. Ebneth, R. Godemann, K. Stamer, S. Illenberger, 
B. Trinczek, E.-M. Mandelkow, and E. Mandelkow, The 
Journal of Cell Biology 143, 777 (1998) 

[8] B. Trinczek, A. Ebneth, E. Mandelkow, and E. Man- 
delkow, Journal of Cell Science 112, 2355 (1999) 

[9] R. Dixit, J. L. Ross, Y. E. Goldman, and E. L. F. 
Holzbaur, Science 319, 1086 (2008) 
[10] B. S. Govindan, M. Gopalakrishnan, and D. Chowdhury, 
EPL (Europhysics Letters) 83, 40006 (2008) 



[11] V. Belitsky and G. M. Schiitz, Journal of Statistical Me- 
chanics: Theory and Experiment 2011, P07007 (2011) 

[12] A. Melbinger, L. Reese, and E. Frey, Phys. Rev. Lett. 
108, 258104 (2012) 

[13] D. Johann, C. Erlenkamper, and K. Kruse, Phys. Rev. 
Lett. 108, 258103 (2012) 

[14] B. Derrida, Physics Reports 301, 65 (1998) 

[15] A. Parmeggiani, T. Franosch, and E. Frey, Phys. Rev. 
Lett. 90, 086601 (2003) 

[16] N. Hirokawa, R. Nitta, and Y. Okada, Nat. Rev. Mol. 
Cell Biol. 10, 877 (2009) 

[17] K. Nishinari, Y. Okada, A. Schadschneider, and 
D. Chowdhury, Phys. Rev. Lett. 95, 118101 (2005) 

[18] Y. Chai, R. Lipowsky, and S. Klumpp, Journal of Statis- 
tical Physics 135, 241 (2009) 

[19] J. Sparacino, P. W. Lamberti, and C. M. Arizmendi, 
Phys. Rev. E 84, 041907 (2011) 

[20] A. Parmeggiani, F. Jiilicher, L. Peliti, and J. Prost, EPL 
(Europhysics Letters) 56, 603 (2001) 



